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ABSTRACT 


The  major  theoretical  efforts  dealing  with  grain-boundary  diffusion 
have  been  reconsidered  with  a  view  to  establishing  their  mathematical  rigor 
and  laying  a  firm  foundation  for  further  work.  The  details  of  Whipple's 
solution  to  the  boundary  value  problem  have  been  reconstructed  and  are  given  in 
this  report.  Fisher's  approximate  solution  and  its  relationship  to  the  exact 
solution  have  been  considered.  A  way  of  evaluating  the  accuracy  of  this 
approximation  has  been  outlined.  The  work  of  Borisov  et  al.,  has  been  shown  to 
be  equivalent  to  that  of  Whipple.  Finally,  an  outline  is  presented  of  compu¬ 
tational  problems  directed  toward  improving  our  understanding  of  the  relationship 
of  grain-boundary  diffusion  to  the  usual  model  and  evaluating  various  approxi¬ 
mations  which  are  frequently  used  in  the  study  of  grain-boundary  diffusion. 
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INTRODUCTION 


Theories  of  grain-boundary  diffusion  invariably  begin  with  the 
formulation  of  a  boundary-value  problem  whose  solution  is  expected  to  describe 
grain-boundary  diffusion.  A  crucial  question  is  then  how  well  does  the  solution 
actually  describe  the  phenomenon.  To  answer  this  two  things  are  needed:  an 
accurate  experimental  study  of  grain-boundary  diffusion  and  an  equally  accurate 
solution  of  the  boundary-value  problem.  Unfortunately  the  exact  solution  of  the 
boundary-value  problem  is  not  usually  possible  in  closed  form  and  in  the  single 
case  where  an  exact  solution  is  known  the  form  of  the  solution  makes  its 
numerical  evaluation  extremely  difficult.  As  a  result  approximate  solutions 
have  been  developed  (in  the  case  being  considered  the  approximate  solution 
antedates  the  exact  solution) .  This  raises  another  question  of  a  purely 
mathematical  nature,  namely  how  accurate  is  the  approximation. 

In  this  report  the  boundary-value  problem  treated  by  Fisher^  and 

(2) 

by  Whipple  has  been  considered  carefully  in  an  attempt  to  resolve  the 
following  questions: 

1.  Is  Whipple's  solution,  in  fact,  exact? 

2.  What  are  the  weaknesses  of  Fisher's  approximate  treatment? 

3.  What  is  the  range  of  validity  of  Fisher' 3  approximation? 

The  first  two  of  these  questions  have  been  answered,  the  second  try- 
pointing  out  that  neglecting  certain  partial  derivatives  is  not  admissible 
under  all  circumstances.  Whipple's  solution  is  shown  to  be  exact  by  means  of  a 
careful  reworking  of  his  analysis.  The  remaining  question  is  partially  answered 
try  outlining  the  calculations  needed  to  justify  Fisher's  solution. 


WHIPPLE *3  SOLUTION 


Statement  of  the  Problem 

The  idealized  situation  studied  by  Whipple^  is  that  in  which  the 
half  space,  y  >  0,  is  filled  with  a  material  of  diffusivity  D  except  for  a  thin 
slab  of  width  2a.  This  slab  is  bounded  by  the  planes  x  *  a  and  x  °  -  a.  The 

'/  ' 

diffusivity  in  this  slab  is  D  (D  >  >  D).  At  time  t  *  0  the  concentration  on  the 
surface  y  *  0  is  suddenly  raised  to  unity  and  maintained  at  unity.  The  problem 
is  to  find  the  concentration  C  ■  C(x,y,t)  elsewhere. 

i 

Let  C  denote  the  value  of  C  in  the  thin  slab, 
equation  says  that 

'  2  '  SC* 

D’c 

in  the  slab.  Outside  the  slab 

D72C  -  .  (2) 

ot» 


Then  the  diffusion 

(1) 


1 


It  is  assumed  that  for  t  >  0,  C  is  a  continuous  function  of  x  and  y. 
This  means  that 


C  -  C  fe.'  x  ■  ±  a. 


(3) 


equation 


9C 

At  the  boundaries  x  ■  ±  a,  ^  is  discontinuous  and  satisfies  the 

d'^L-d^  . 

3x  3x 


(h) 


The  derivatives  involved  in  Eq.  (h)  are  one  sided  derivatives  and  the  prime  (or 
lack  of  one)  indicates  the  direction  in  which  they  are  taken. 

On  the  surface  y  -  0, 


When  t  a  0, 


C(x,0,t)  ■  1  when  t  >  0. 

(5) 

C(x,y,0)  -  0  if  y  >  0. 

(6) 

Approximate  Boundary  Condition  at  x  8  a 

An  approximate  boundary  condition  on  the  surface  x  a  a  vail  now  be 

derived. 

Because  of  the  symmetry  of  the  situation  C*  is  an  even  function  of  x. 
Hence  it  can  be  written  as  a  power  series  in  x  with  the  odd  powers  of  x  missing. 
The  coefficients  are  functions  of  y  and  t. 

2 

c'(x,y,t)  -  C0(y,t)  +  y  C2(y,t)  +  ...  .  (7) 

When  the  derivatives 


y 


+  ...  ,  and 


y 


are  substituted  in  Eq.  (1)  and  terms  involving  second  or  higher  powers  of  x  are 
dropped,  one  obtains 

i 

D 


30o  •' 

~T  *  C2I 
ay 


3C( 

3t" 


(8) 


3t 

2  ' 
3_C_ 

ay2 

aV 


fo 

3t  + 

'V 


C2  + 


2 


for  points  in  the  slab  (in  the  applications  a  is  a  very  small  number). 

Again  dropping  powers  of  x  which  are  second  or  higher  there  results 

i 

C  ■  Cq  and 

D  fj  ■  o'  aCg  at  x  -  a, 

because  of  Eq.  (3)  and  Eq.  (U)  and  the  fact  that 

ir "  +  •••  • 

1 

Now  since  C  ■  CQ  at  x  ■  a,  and  the  partial  derivatives  with  respect 

to  y  and  t  on  the  surface  x  ■  a  depend  only  on  the  common  value  of  C  and  C„  on 
this  surface, 

3C  .  £|o 

at  at 


2  2 

ay^  ay^ 


at  x  ■  a. 


Substituting  all  of  these  quantities  in  Eq.  (8)  gives 

D'  H + « I  ■  § for  *  ■ a-  (u) 

ay 

The  differential  Eq.  (2)  can  be  used  to  write  this  last  equation  in  a 
more  usable  form.  Putting  x  ■  a  in 

ax2  ay2  at 

which  holds  for  all  points  outside  the  slab,  and  combining  this  with  Eq.  (U)  gives 


«'  R&ml]L  n]  ac 

,  2  "  a  ax  D  at 


for  x  ■  a. 


This  is  the  form  in  which  the  approximate  boundary  condition  will  be  used. 


The  Fourier-Laplace  Transform 


In  order  to  solve  the  diffusion  equation  the  function  C(x,y,t)  is 
transformed  by  both  a  Fourier  and  a  Laplace  transform.  The  function  defined  is 


Kx,H,X) 


exp(-\t)dt  [  sin(|iy)  C(x,y,t)dy. 


3 


For  purposes  of  discussion  the  intermediate  function 


where  a  is  the  real  part  of  X,  and 


lim  2 W  exp(-at)  exp  &J  -  0. 
t-^«»  _  / 

This  also  shows  that  4(x,n,X)  has  a  bound  which  is  independent  of  x  and  p,  (but 
not  X). 

When  the  function  \Kx,p,X)  is  obtained,  it  will  be  seen  to  possess  the 
correct  properties  so  that  the  inversion  formula 

Y  +  iR 

6(x,p,t)  "  2^  lim  f  exp(Xt)  vKx,n,X)dX 
Rf*#0  \  -  ill 


holds  for  y  >  0.  From  the  properties  on  C(x,y,t)  it  will  be  shown  that 


00 

C(x,y,t)  -  |  f  6(x,p,t)  sin(py)dp. 


Transformation  of  the  Differential  Equation 


When  C(x,y,t)  is  transformed  by  the  Fourier-Laplace  transform,  the 
o  o  7  7 

derivatives  8C/8x  and  8  C/8x  are  transformed  into  di'/dx  and  8 i/dx  ,  respectively. 

2  2 

The  derivative  3  C/dy  will  now  be  transformed.  First  it  must  be 
observed  that 

lim  C(x,y,t)  ■  0,  and 
y-+oo 

lim  |£  -  0. 

y-  * 


Then  the  integration 


r 

sin(py)  ~ 


dy  can  be  performed  to  give 


5 


«  /  2  \  \  r  \ 

f  sin(py)  p4  dy  ■  sin(py)  "  H  /  c°s(py)  dy 

1  3y  /  *  „  .  n  Jn  >' ** 


y  -  0  0 


2  ( 

-  p  cos((jy)  C(x,y,t)  -  | x  I  sin(py)  C(x,y,t)dy 

-  y  -  0  J0 


c,0,t)  -  \?  ^ 


C(x,y,t)  sin(py)dy. 


Observing  that  C(x,0,t)  *  1,  multiplying  this  last  expression  by  exp(-Xt),  and 
integrating  gives 


exp(-\t)dt  sin(py)  — -  dy 

-in  M 


“  CO  CD 

“  p.  I  exp(-Xt)dt  -  p2  f  exp(-Xt)dt  f  C(x,y,t)  sin(py)dy 

•^A  J  A  “'A 


£  -  p2  '[(xjij.jX)  . 


Thus  — -  is  transformed  into 

ay 

ix  2 

-  -  >Kx,n,x)  . 

9° 

To  transform  the  integral  exp(-Xt)  ^  dt  is  evaluated.  Because 
C(x,y,0)  =  0, 

®  I  \  I  t  ■»  «, 

f  exp(-Xt)  ~  dt  3  exp(-Xt)  C(x,y,t)  +  X  f  exp(-Xt)  C(x,y,t)dt 

J o  J  t  «  0  0 


© 

■  X  I  exp(-Xt)  C(x,y,t)dt. 

■'a 


Multiplying  by  sin(py)  and  integrating  gives 


j  sin(py)dy  J  exp(-Xt)  |-  dt  ■  X  f  sin(py)dy  j  exp(-Xt)  C(x,y,t)dt. 

*  **/N  •#  A 


6 


Thus  — *  is  transformed  into 
ot 


X>K  x>M>>X)  • 


From  these  results  it  follows  that  the  differential  Eq.  (2)  is 
transformed  into  the  equation 


Rewriting  this  gives 


D  ^  -  |i2i  B  Xt. 

3x 


4- 


The  boundary  condition  Eq.  (12)  transforms  into 

D'  4  - » i "  (r  ■  ^  M' tor  *  ■ a- 

ax  / 


Solution  of  the  Transformed  Equation 

The  function  'Kx,n,X)  is  to  be  regarded  as  a  complex  function  of  the 
real  variable  x  for  fixed  p  and  X.  Then  Eq.  (16)  is  essentially  an  ordinary 
differential  equation  in  x  with  complex  coefficients. 

The  solution  of  Eq.  (16)  is 

>Kx,p,X)  -  A  exp  xVp2  +  +  B  exp  -  xVp2  +  ^  +  E  ,  (3 

»  / 

where  A,  B,  and  E  are  functions  of  X  and  p  and  the  radical  indicates  that  the 
square  root  with  positive  real  part  be  taken.  Since  t  is  a  bounded  function  of 
x,  we  can  reject  the  term 

A  exp  x\/  p^  +  —  . 

This  leaves  B  and  E  yet  to  be  determined.  Because  the  term 

B  exp  -  xVp2  +  “| 

is  a  solution  of  the  homogeneous  equation 

ax  / 


it  can  be  concluded  that 


-  i***5 


7 


and  therefore 


2  .  \\  ‘ 

D 


(19) 


The  function  B  must  now  be  adjusted  so  that  satisfies  the  boundary 
condition  Eq.  (17)  for  x  ■  a.  Substituting  % 


2! 

0X 


*  -  B 


V,2  +  I  exp  -  xV^I 


and 


^|-B 

dx^ 

into  Eq.  (17)  when  x  B  a  gives 
B 


2  x  x 

a  +5 


exp 


-\[~2  X' 
“  xV  p,  +  - 


" 

/  \ 

H 

I  ^  \ 

—  1 — 1 

1 

"ala 

1 

D 

1  2  V 
i*  +5, 

A  D,/ 2  X 

+  aV^  +  5‘ 

X 

exp 

X|a 

+ 

1 

u  | 

l  ) 

- 

\  I 

l  / 

Xu 


\  2  x 
x>*  +i 


Solving  for  B  one  obtains 


D 


D 


-  -  1 


— \  /  p  \ 

exp  ayii  ♦  ^ 


2  ,  Dn/  2  .  X  . 

H  +  i  +  d  + 


Therefore  \Kx,[i,X)  is  determined  uniquely  to  be  ^  +  tg  where 


1  ^  2  A  x\ 
xd  ♦  d 


and 


I  » 

JL.  i 

Id  / 

H  exp 

-  (x-a)^2  ♦  i 

2  *  x^ 

+D 

D'  *  +  +  5  +  x 

Inversion  of  t. 


(20) 


0^  will  be  used  to  denote  the  inverse  of  The  formula  for  is 


y  ♦  iR 

4  sin(py)dn 

lin  f  oxplXAL  dX 

/ 

0 

•iU-  -V  -  iR  »*  +  d)  - 

where  y  >  0*  The  value  of  the  quantity  in  the  brackets  is  easily  seen  to  be 
independent  of  y  as  long  as  y  >  0.  This  will  be  discussed  later. 

8 


When  the  substitutions} 

y  ■  T)VSt ,  \i  ■  “7=  and  X  ■  7-  } 
VDt  t 

are  made,  reduces  to 

°i  ■  7:  / 

n  x  Jn 


sin(jj,  rj)dp. 


y  +  iR 

t 

1  " 

lim  f 

.  ^  -y  -  iR 

exp 

X 

dX 

•2 

V  X*) 

1 

X 

Since  the  value  of  was  independent  of  y  before  the  substitution 

the  last  form  is  also  independent  of  the  choice  of  y  as  long  as  y?  0.  The 
primes  can  be  dropped  and  C1  can  be  written 


00 

r 

■ 

y  +  iR 

- 

j  p,  sin(piT])dp. 
0 

lim 

R-*aO 

m 

f  exp  X 

dX 

Y  -  IR 

X 

For  the  purposes  of  evaluating  the  contour  integral,  X  will  be 
written  as  x  +  iy.  (There  should  be  no  confusion  of  this  x  and  y  with  the 
original  x  and  y  of  the  problem.) 


The  integrand, 


exo  X 


is  analytic  for  all  values  of  X  except  for  X  =  0  and  X 
simple  poles.  Therefore  the  integral  of  this  function 
which  encloses  both  poles  has  the  value 


cxp(-(i2) 


M- 


2 

3  -  (i  where  it  has 
about  a  closed  contour 


since  the  residue  at  X 


0  is  ^  and  the  residue  at  X 


Vi.-sabil 

M- 


Let  the  closed  contour  to  be  used  consist  of  the-  line  x  3  y  between 

|  IQ 

the  points  y  -  iR  and  the  point  y  +  iR  and  the  part  of  the  circle,  X  3  R  c  , 

•  2 

Vr  +  y  .  Call  this  contour, 
1  1  2 

Suppose  R  >  p.  .  Then 


which  is  to  the  left  of  the  line  x  3  y,  where  R 

1 

as  shown  in  Figure  1,  I*.  Note  that  as  R-t*  so  does  R 


I 


exo  X 

^  *  N 


dX 

X 


2ni 


2 

1  -  exp(-u  ) 
2 


9 


10 


The  length  of  the  curve  P^P.^  is  less  than  2nR  j  therefore 


/ 

_|=ai_dX 

p  p  p 

(p^  +  \)\ 

2n  exp  y 
(R*  -  p2) 


This  integral  can  be  seen  to  approach  zero  as  .  Hence 

Y  +  iR  r  o 

r  _ 1  j-v  N  .  - ( 


r  +  iR  r  2 

lia  {  4SLL.  *  o  2ni  Lz-SfhLL 

ft*.  {  -  ili  (|1  +x>  n 


This  permits  to  be  written: 


\a  n  f  sinM  1  ‘  exP("l^2)  ^  *  f 

o  L  J 

f  dp.  -  -  f  dp  f  cos(|xr]  )  exp(-p2)dp 

u  -'0  JQ 

T] 

,  f  *2  , 

■  1  -  *^=  J  gM-  \“)dT)  ■  1  -  erf  j)  •  erfc  |  . 


cos(  [XT]  )  exp(-p  )dp 


These  last  integrals  r.ay  be  evaluated  from  a  standard  table  of  integrals, 
for  example,  formulas  lj06  and  la32  from  the  table  of  integrals  of  the  Handbook  of 

Chemistry  and  Physics^ 


Inversion  of  i 


Let  denote  the  inverse  of  >!-  Then 

•  r  V  +  iH 

C  =  /  sin(py)dp  lira  /  exp(\t)t>  dX  * 

nhJ0  MY-1R  J 

where  y  >  0.  It  will  be  seen  later  that  C?  is  independent  of  the  choice  of  y 
as  long  as  y  >0. 


The  following  substitutions  will  now  be  made: 

K  "  ~  »  *1  "  >  a  "  "TTi  A  =  ~  ,  p  «  ~  *  U-l)a,  p  c  and  X  -  . 

^Dt  -fDt  m  u  -fDt  ^Dt 


11 


i8  reduced  by  these  to 


tYDt  (A-l)p*  exp[  -  ^l/p  + 


/  »2  »w.  *2  1,1  *2  »  ». 
(H  +  X  )(4i  £  X  +  X  , 


When  this  is  put  into  the  expression  for  C.  and  the  substitution 
completed  one  obtains 


Y  ♦  iR  f~o - 

c  .M  Lainl^  lin  f  12^— 

o  L^“  -V  -  in  (n2  ♦  x)(*i2 1  i^n2  ♦  x  +  x). 


after  dropping  the  primes.  Note  that  y  >  0  and  this  expression  for  C  inherits 
the  property  of  being  independent  of  y»  ^ 


It  is  desirable,  however,  to  take  the  contour  integral  along  the  line 

2 

x  ■  -  p  .  (Again  X  is  written  x  +  iy  for  the  contour  integration.)  The  point 
2 

X  ■  -  p  is  a  branch  point  of  the  integral  and  must  therefore  be  avoided.  The 
contour  which  will  be  used  is  shown  in  Figure  2  and  is  denoted  by  T(R). 


i 

I 

j 

I 

1 


The  curved  part  of  r(R)  is  a  semi-circle  of  unit  radius  with  center 

P  P 

at  the  point  -  jjt  ,  while  the  rest  of  r(R)  consists  of  the  line  x  ■  -  \r  between 

o  2 

the  points  -  p  -  iR  and  -  p  +  iR. 


The  function 


V77 


X  is  analytic  for  all  values  of  X  where  X  is  not 


2  2 

real  and  all  real  values  of  X  greater  than  -  p  .  The  only  zero  of  (p  +  X)  is 

2  ?  i  /  p 

at  X  ■  -  ja  .  The  function  Apr  +  -  yp  +  X+  Xis  analytic  at  all  points  where 

V2  a 

p  +  X  is.  It  cannot  be  zero  if  X  has  an  imaginary  component  because  a  >  0 

and  Im(V  p^  +  X)  has  the  same  sign  as  Im(X).  It  cannot  be  zero  for  real  X 
2 

greater  than  -  p  because  A 
case. 


>  1  and  W  +  X  is  a  positive  real  number  in  this 


Therefore  the  function 


exp(X  -  rin'  +  X) 


(p2  +  X)(Ap2  +  -  ViJ+X  +  X 


f(X) 


is  analytic  at  all  points  cn  T(R)  and  all  points  to  the  right  of  f(R).  From 
this  it  follows  that 


/ 

r(R) 


Y  +  iR 

f(x)dX-  f  f(X)dX 

Y  ~  iR 


f(X)dX  + 


f(X)dX  . 


V3 


P7P1 


Along  the  segment  P^,  *  *  *  +  iR  anc! 

|p4  +  X|  >  d, 

,  !a  p2  + 1  V^2  +  x  +  xl>R 

(because  Ir.( V p2  +  X)  has  the  same  cvr.  as  Ir(X)),  and 


I c ::r  X  -  ?->Jp2+  X  |  <  exp  y 
(because  Icc(V  p2  +  X)  >  j).  This  gives  that 


f(X)dX 


F  P 

2*3 


ex?  Y 

r,2 


(Y  +  |i  ) 


and  this  goes  to  zero  as  R->«»  .  A  similar  argument  holds  for  PyP^*  These 


13 


statements  show  that 


^  I 


f(\)d\ 


Y  +  1R 


Y  -  iR 


f(X)d\  . 


The  transformation  v  ■  V|i^  +  X  will  now  be  performed  on  the  integral 
along  f(R),  where  again  the  radical  indicates  the  square  root  with  positive  real 

part.  This  transformation  takes  T(R)  into  the  contour  f  (R)  shown  in  Figure  3. 


This  contour  consists  of  parts  of  the  lines  y  ■  x  and  y 
fourth  of  the  unit  circle  with  center  at  the  origin.  Thus 


and  one- 


J  exp(X  -  Hi*  +  X)dX _  _  2  exp(-  /)  f  exp(/  -  vg)dv 

r(R)  +  +  “  Vp-2  +  *  +  x)  (A  -  1)  (|12  +  i  v  +  — 

P  A  -  i 

where  p  -  (A  -  l)a  as  before. 


2  exp(- 


2  .  1 


lh 


Since  the  integrand 


is  analytic  for  all  v  in  the  half  plane  x  >  0,  r*(R)  can  be  replaced  by  the 

t  I  P 

contour  f  (S)  shorn  in  Figure  U,  where  S  ■  . 


This  gives  the  following  expression  for  C ^ 


00 

|-  J p.  exp(-n^)  sin(pr))dp, 
1  0 
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Along  the  contour  r  (S)  the  following  substitutions  are  valid: 
M-  sin(nri)  ■  -  cos(p-)),  ^2! — ^  ex p(- 


(because  Re(v)  >  0  on  T  (S)  ),  and 


X  2  v  v  p  \i  o 

- —  -exp^+T  +  pY  exp  -  H  +  o  +  pi  a  d a  . 

+  '  '  1  L.'  'J 
2  v  t i 

(because  p,  +  |  has  positive  real  part  along  T  (S)  ). 

Applying  these  substitutions  to  gives 


CO 

2  8  f 

-  -r-  -r-  I  cos(p.ri)dp,  lim 

n  i  011  <  S->o 


where 


,  CO  CO  r 

Q  ■  J  d v  J  exp(-p,2a)da  J  exp  pj  “  v(£  + 


r  (S) 


The  limit  of  Q  as  $+*>  is  actually  an  improper  real  integral. 
Therefore  it  is  permissible  to  interchange  this  limit  process  with  the  various 
integrations  to  get 

«  go  co  r  f  . 

°2  *  ”  "It  8tq  J  ^  J  da  /C0SW  exp(-^2a)d|i  lim  J  exp  v2 

n  1  ?  io  Ls**r,,(s) 


The  last  integration  must  be  divided  into  two  cases.  The  case  where 
o  >  A  gives  an  integral  of  the  form 


lim 

S+ao 


Gv)dv 


where  F  and  G  are  positive  real  numbers.  The  integrand  is  analytic  for  all 

ii.  . 

values  of  v.  Therefore  the  integral  around  a  closed  path  is  zero.  If  r  (S) 
is  closed  by  the  vertical  line  from  S  +  iS  to  S  -  iS  one  obtains  the  contour  shovm 
in  Figure  $. 


values  of  v 


Along  the  line  segment  P^P^,  v  ■  S  +  iy  where  -  S  <  y  <  S.  For  these 


|exp(-Fv2  -  Gv)|  *  exp 


-  F(S2  -  y2)  -  03 


<  exp(-GS). 


This  drives 


/ 

Vi 


2 

exp(  -  Fv  -  Gv)dv 


<  2S  exp(-GS), 


and  this  goes  to  zero  as  .  Thus 


lim 
S-+  ap 


f 


r"(s) 


exp( -Fv 


2 


Gv)dv  s  0. 
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The  case  where  o  <  t  gives  an  integral  of  the  form 


lim  I  exp(Jv^  -  Kv)dv 
r"(s) 


where  J  and  K  are  positive  real  numbers.  This  will  be  shown  to  be  equal  to 

Y  +  iS 

lim  f  exp(Jv^  -  Kv)dv  . 
s*“  ■y  -  iS 


To  do  this  it  must  be  observed  that 


1 


exp(Jv  -  Kv)dv 


ru(s) 


differs  from 


Y  +  iS 

f  exp(Jv^  -  Kv)dv 

^Y  -  iS 


by  the  contribution  from  the  horizontal  segments  and  shown  in  Figure  6 


* 


For  the  integral  along  the  segment  Pj^,  v  ■  x  +  iS  where  Y  <  x  <  s* 
Vherl  v  is  on  the  segment  of  this  line  where  Y  <  x  <  S/2> 

s2 

|exp(Jv2  -  Kv)|  -  exp[j(x2  -  S2)  -  Kx]  <  expjj 


a  -  S2 


exp 


-¥ 


When 


v  is  on  the  segment  of  where  S/2  <  x  <  S, 


|exp(Jv2  -  Kv) |  ■  exp 


(j(x2  -  S2)  -  Kx 


1  KS 

c  <  exp 
J 

r  2 

Therefore 


I  exp(Jv2  ”  Kv)dv 

<  s 

exp 

-  ^iJL  |  +  exp  - 

4 

PUP5 

and  this  goes  to  zero  as  S+*>  . 

This  now  leaves  for  Cg  the  expression 


«  A  90 

.  JL.  SL.  f  tig*  f  do  I  cos((xt])  exp(-  ^2e)Qdp, 
n2i  ^  |  0 


where 


Q  =  lim 

Y  +  iS 
( 

exp 

h-a  2 

Pi v  -u(5  T 

S-+» 

y-  is 

- 

dv  . 


let  V  be  written  r  ♦  i y  lor  the  evaluation  of  Q.  Then  dv  -  idy  and 
Y  +  iS 


exp(Av  ~  Bv)dv 


;Y  -  iS 


3 

. 

J  6Xp 

-S 

a(y2  -  y2  +  2yyi)  -  b(y  +  iy) 

i  exp(AY2  -  By)  j  exp(-Ay2)  cos 

-S 


dy 


(2Ay  -  B)y 


<37 


(21) 
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i 


This  gives  the  final  solution 


c  -  01  +  c2 


where 


C1  =  erfc  j 


Discussion  of  Inversion  Formulas 


In  the  half  plane  2c (\)  >  0  both  t,  - - - - and 

|j  -  3)l*  exp[-  Vt*2  +  3  U'a)] 

|F * 5)  IF  F  ♦  +  x) 

are  of  order  l/X2.  When  the  real  part  of  X  is  not  negative,  )pp.2  +  x|  >  |x| 
and  It.  X2|  <  d||xI  .  Similarly  when  Re(X)  >  0,  I D *  p.2  +  +  7-  +  X|  >  |x| 

1-n  I  if*1  '  T*  "  a  j  — 


because 


LD  !~2 - T 

rvn  +p 


has  its  real  and  imaginary  parts  the  same  sign  as  those  of  X. 


!t2x2l  <  [§-  - 1  d|hI. 

u  1 


When  Re(X)  >  0  both  t.  and  t«  are  analytic  functions  of  X  and  are  both 

12  f  o 

real  if  X  is  real.  Therefore  by  a  theorem  in  Churchill^  '  the  inverse  of 
t^  3  t^  +  t2  is  given  by  the  formula 


Y  +  iR 


l  f 

x,n,t)  a  lim  t  exr 

J  -n 


exp(Xt)  dX, 


Y  -  iR 


for  y  >  0>  and  is  independent  of  y- 
It  was  shown  before  that 


J  C(x, 


y,t)dy  exists.  Therefore  since 


C(x,y,t)  is  assumed  to  be  a  continuous  differentiable  function  of  y,  one  can 
write ^ 


«  -  r* 

2  f  1  1  1  t  2  r 

C(x,y,t)  -  -  J  sin(py J  G(x,y  ,t)sin(|jy  )dy  -  -  J  sin(py)a(x,ji,t)dji. 
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Let  C(x,y,t)  be  a  solution.  Then  0(x,p.,t)  and  'Kx^X),  determined  by 
C(x,y,t),  must  exist.  >Kx,[i,\)  is  then  uniquely  determined  and  uniquely  inverted 
to  0(x,u,t).  But  0(x,|i,t)  has  a  unique  inverse  also  hence  there  is  only  one 
possible  choice  for  C(x,y,t). 


FISHERY  SOLUTION 


Statement  of  the  Problem 

(2) 

The  idealized  situation  studied  by  Fisher'  '  is  the  same  as  that  of 
Uhipple.  The  half  space,  y  >  0,  is  filled  with  a  material  of  diffusivity  D 
except  for  a  thin  slab  of  width  5.  This  slab  is  bounded  by  the  planes  x  a  0 

I.  i 

and  x  e  -  6.  The  diffusivity  in  this  slab  is  D  (D  >  >  D).  At  time  t  *  0 
the  concentration  on  the  surface  y  =  0  is  suddenly  raised  to  unity  and  maintained 
at  unity.  The  problem  is  to  find  C  ■  C(x,y,t)  elsewhere.  Fisher  refers  to  the 
plane  y  0  as  the  free  surface. 

The  diffusion  equation  holds  for  the  concentration  G  outside  the  thin 
slab.  That  is 

d72c  =  h,  (26) 

for  C  outside  the  thin  slab. 


The  Differential  Equation  in  the  Slab 

Let  the  thickness  of  the  slab,  6,  be  very  small  so  that  variations  in 
C  across  it  are  negligible.  Consider  a  rectangular  parallelopiped  in  the  slab  of 
length  Az  in  the  z-direction,  of  'width  Ay  in  the  y-direction  and  of  width  slightly 
greater  than  6  in  the  x-direction.  Let  this  be  oriented  in  the  slab  with  two  of 
its  faces  .-just  outside  the  slab  and  parallel  to  the  slab.  Let  P  be  the  point  at 
the  center  of  this  parallelopiped.  Let  dC/dt,  denote  the  value  of  that  derivative 
at  the  point  P.  Then  the  flow  of  material  into  this  volume  is  approximately 
6( Ay )(Az)  dC/dt. 

This  can  also  be  written  as  the  sum  of  the  flow  through  the  sides  of 
the  parallelopiped.  In  the  z-direction  there  is  no  flow.  In  the  x-direction 
there  is  the  flow  F  just  outside  the  surface  x  B  0  (and  its  negative  just 

outside  the  surface  x  ■  -  5).  Therefore  the  amount  of  material  coming  into  the 
volume  in  the  xrdirection  is 


-2(Fx)q  (Ay)  (Az), 


where  the  subscript  zero  indicates  that  F  is  evaluated  just  outside  the  slab. 
Of  course  by  Pick’s  law  (outside  the  slab) 
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The  flow  at  the  point  P  in  the  y- direction  is  F  and  so  the  amount 

y 

of  material  coming  into  the  volume  in  the  y-direction  is 


F  .fik 

7  ay  2 


Uz)6  - 


ap  A 
F  +  —X  k 
y  ay  2 


3F 

(Az)5  ■  -  ^  (Ay)U z)6. 


t 

Noting  that  inside  the  slab  F  ■  -  D  ■—  and  equating  the  expressions 

y  oy 

for  the  amount  of  material  flowing  into  the  volume  gives 

-2(Fx)0  ( Ay)( Az)  -  ^  (Ay)(Az)6  «  |£  (4y)(Az)6 

or 


2D  /6C\  .  '  a2c  =  SC 

"5  tel  *D  75  at 

^  oy 


(27) 


for  the  concentration  in  the  slab. 


Fisher »s  Approximate  Solution 

At  this  stage  the  problem  is  identical  to  that  solved  by  Whipple. 
Equations  (26)  and  (27)  are  just  the  same  as  Eqs.  (2)  and  (11). 

Equation  (26)  is  simplified  by  the  observation  that  the  concentration 
of  the  diffusing  material  will  be  much  greater  in  the  slab  than  outside  because 
of  the  much  higher  diffusivity  inside.  This  means  that  near  the  surface  x  “  0 
and  far  from  the  surface  y  ■  0  the  flow  of  material  is  normal  to  the  surface 
X  ■  0. 

In  terms  of  the  parameters  of  the  problem  this  says  that 


From  this  it  is  concluded  that 


Equation  (26)  is  therefore  reduced  to 


D 


dC 

at 


for  points  outside  the  slab. 


(28) 


23 


Let  the  following  substitutions  be  made: 


Dt  x  y 

Equation  (28),  written  in  terms  of  these  dimensionless  parameters,  becomes 


a  c  n  sc 

3^  3tl 


(29) 


for  points  outside  the  slab  and  Eq.  (27)  becomes 


2|SCJ 

*  2  at1 
^  i 


'9xi  o 


(30) 


for  points  in  the  slab. 


Fisher  solves  these  two  equations  by  a  numerical  method  and  uses  the 
numerical  solution  to  draw'sone  conclusions  about  C(x,y,t).  He  observes  that  C 
rises  in  the  slab  at  a  rapidly  decreasing  rate.  Because  of  this  the  value  of  C 
outside  the  slab  can  be  treated  as  though  the  value  of  C  in  the  slab  is 
independent  of  the  time.  Let  <p(y^,  O  be  the  concentration  in  the  slab. 

The  assumption  that  <p(y  ,  t, )  is  independent  of  the  time  gives  for  a 
solution  of  Eq.  (25)  1 


s J 


p(y1» 


* 

i 


and  this  takes  on  the  value  <*>(y^,  t^)  when  x^  =  0. 


(31) 


At  this  point  Fisher's  reasoning  is  difficult  to  follow.  He  says  that 


yr  in  the  slab 
ot 


SC, 

>5,  at 


outside  the  slab  for  most  of  the  time.  From  this 


max 


he  concludes  that  the  value  of  p(y, ,  t, )  can  be  derived  approximately  from  Eq.  (30) 

SC  ^  ^ 

by  assuming  B  0. 


This  assumption  applied  to  Eq.  (30)  is  that 


However 


1 

Vn"l“ 


t^)  exp 


(32) 


2ii 


t 

I 


for  points  outside  the  slab,  therefore 


|f~)  -  -  —=  tJ  • 


Putting  this  last  egression  into  Eq.  (32)  gives 

V  "  3y‘ 

which  is  essentially  an  ordinary  differential  equation. 

The  solution  to  Eq.  (33)  which  is  bounded  in  y  and  satisfies  the 
condition  that  0  ■  1  when  y^  ■  0  is  1 


(p(y1,  tx)  -  exp 


V2  ' 
7K  tdA 


This  gives  the  final  formula  for  C. 

V2"  7l 


G  =  exp 


M ♦  iA 


erfc 


*1 


Tjt1 


or  in  terns  of  x,  y,  and  t 

C  *  exp 


V7: 


Vs’1'  jjj-  (nDt)1^ 


erfc 


x 


2^t 


(3k) 


DISCUSSION  OF  VALIDITY  OK  FQfttl JUS 

The  solution  derived  by  Whipple,  Eq.  (25),  is  the  exact  solution  to  the 
diffusion  equation  with  the  approximate  boundary  condition  Eq.  (11),  the  boundary 

conditions  Eqs.  (k),  (5),  (6),  and  the  conditions  that  ~  and  C  both  go  to  zero  as 

y-+  00  •  His  solution  has  the  disadvantage  of  being  a  numerical  integration. 

Fishers'  formula,  Eq.  (3k),  on  the  other  hand,  is  ideal  for  computation 
if  it  is  correct.  The  reasoning  of  Fisher  will  now  be  discussed. 

First  the  simplification  of  the  diffusion  equation  to  the  form 

§-d4  (35) 

ox 
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would  have  to  be  verified.  This  can  be  checked  by  using  the  exact  solution  of 

a  p  5  0  p  p 

Whipple  and  comparing  31] /ax  and  3  C/3y  .  If,  as  seems  not  unlikely,  3  C/3y 

2  2 

is  negligible  compared  to  3  C/3x  for  values  of  x  and  y  of  interest,  then  the 
concentration  in  that  region  would  satisfy  Eq.  (35). 

The  next  step  is  difficult  to  check.  Fisher  concludes  that  ^(y^jt^) 
is  a  slowly  varying  function  of  time  and  therefore 

„  _  ,  j.  \  l  *1  1 


C  -  <p(ylf  tx)  erfc 


2Vt "I 


(where  cp  is  still  an  unknown  function)  satisfies  the  equation 


This  is  true  if 


3C  „  3  G 


2L  erfo  Jl' 

*1  2^ 


is  small  compared  to 


2V3T  t 


*ll 

cp  exp  -  Tjq  • 


The  only  way  to  verify  this  seems  to  be  to  look  at  Fisher’s  final  solution  and 
determine  if  this  is  so. 


One  observation  seems  pertinent  here.  If  the  diffusion  equation  can 
be  reduced  to  Eq.  (35)  then  this  last  step  is  a  necessary  condition  for  Fisher’s 
solution  to  be  a  solution  to  the  diffusion  differential  Eq.  (26).  If,  on  the 
other  hand,  Eq.  (35)  is  not  valid,  then  Fisher’s  solution  might  still  be  a  good 
approximate  solution  even  if  this  last  test  fails. 

The  next  step  is  again  one  which  can  be  checked  by  means  of  Whipple’s 
solution.  Fisher  uses  the  equation  in  the  slab. 


2i§y 

dXl'  0 


to  determine  cp(y^,  t^).  This  requires  that  3C/3t^  in  the  slab  be  small  compared 


to  2(3 C/ax^)0.  Since  wisher’s  solution  satisfies  Eq.  (36)  exactly,  this  last 
condition  is  a  necessary  one  for  it  to  satisfy  the  boundary  condition  Eq.  (27). 


DISCUSSION  OF  BORISOV,  GOLIKOV  and  ljubovs  work 

The  following  is  a  partial  study  of  the  derivation  of  Fisher1 s  solution 
from  that  of  Whipple's  by  Borisov  et  al.^) 

Consider  the  formula  given  by  Whipple  for  the  concentration  C  in  the 
grain  volume: 


c  -  c1  +  c2 


where 


'1  n 


1  -  exp(-p  ) 


sin(|XT}) 


_c|x 


and 


2  2  2 
c2  *  ”2“  /  |i  exp(-n  )  sin(|i, 

n  i  < 


lim 
S->  •  j 


f  -  cxp(v^  -  v£)dv 


ft,  V 

r  (s) 


2  V  V 

*  +P  +  FI 


V  J 


y 

It  should  be  noted  that  in  the  convention  used  by  the  Bovisov  et  al.  t)  ■  — *  , 

V~Dt 

as  in  Y/hipple's  paper,  but  £  ■  — ~  where  x  is  the  distance  measured  from  the 

■^Dt 

grain  boundary.  (Whipple  has  £  =  .)  All  other  symbols  in  C.  and  C_  are 

^Dt 

the  same  as  in  Whipple 1 s  paper. 

The  above  expression  is,  of  course,  not  the  final  form  of  Whipple's 

section. 


It  is  desired  to  evaluate  the  contour  integral  of  C„  along  the  y-axis 
.  ",  ,  ,  .  d 
instead  of  T  (S).  More  precisely,  let  F(5,e;  be  the  contour  which  consists  of 

the  part  of  the  y-axis  between  the  points  -  i5  and  iS  excluding  the  segment 
between  -  ie  and  ie.  (It  is  assumed  that  S  is  a  very  large  positive  real  number 
while  e  is  a  very  small  positive  real  number.)  The  remaining  part  of  f(S,e) 

iG 

will  be  the  half  of  the  circle  v  *  ee  which  is  to  the  right  of  the  y-axis. 

This  contour  is  shown  in  Figure  7. 
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-iS 


FIGURE  7. 


It  should  be  noted  that  the  x  and  y  of  Figure  7  and  the  y  to  be  used 
as  a  variable  of  integration  on  contour  are  not  the  original  x  and  v  of  the 
problem.  J 


Let  f(v) 


exp(v2  -  v£) 


A 


4* 


V 

A-l 


The  function  f(v)  is  analytic  on 


T(S,e)  and  at  all  points  to  the  right  of  T(S,c).  This  is  because  the  zeros  of 
2  v  y*- 

**  +  1  +  A-l  are  a11  40  the  left  of  the  y-a3ds  except  for  the  possibility  v  ■  0. 
f(v)  is  bounded  by  exp(v2  ■  v?)  as  .  Therefore  by  an  argument  very  much 
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■M 


like  that  used  .just  before  Eq.  (21)  it  can  be  concluded  that 


lira 

n 

0-+m 


i 


exp(v  -  v£)dv 


exp(v  -  vg)dv 


<M 

■> 

■> 

CM 

—  xxiu  r 

s+«  J 

2  A  v  .  v2 

r 

V 

11 

r(s,e) 

Therefore  it  is  required  to  evaluate 

exp(v  -  vC)dv 

P(S,e) 


I 


2  v  v‘ 

^  +  p +  Vi 


This  integral  is  clearly  independent  of  e  and  so  it  is  permissible  to 
find  the  limit  as  . 

First  on  the  section  of  r(S,e)  which  consists  of  the  half  circle 
v  =  ee^,  -  -  <  8  <  ^  ,  the  integral  becomes 


r 

a 


n/2 


exp(e2  e210  -  £ee10)d0 


-n/2 


2  eei0  e2e2i0 

“  +nr  +  ^r“ 


As  e->  O  this  integrand  approaches  l/|i  uniformly  in  0,  therefore  the  integral 
approaches 


1 T 

-n/2 


®  as  e-*.0  . 


For  the  section  of  r(S,e)  which  lies  along  the  y-axis  the  substitution 
v  *  iy  is  made.  On  this  section  the  integral  becomes 

f  exp(-y2  -£yi)  dy  f  exp(-y2  -  £yi)  dy 


2  y2! 

si 

v,l 

+ 

_  » 

!  2  y2 

^  -j=i 

1 

“PiJ 

-  S 

If  y  is  replaced  by  ”  y  in  the  first  integral  one  obtains 
S 


iy  y 
+  P 


/. 


e*p(-y2  -  gyi)  exp(-y2  +  ?yi) 


m2  .  jdl 

+  1L.  1 

C\J_ 

.JJL 

y 

*  A'1 

P  1 

^  A-lj 

P 

dy 
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-  2i  exp(-y  ) 
J 


-  H2  sin(gy)  -  |  cos  gy 


iM 


dy 

y 


In  this  form  it  is  clear  that  there  is  a  limit  as  e  -*0  and  this  limit  is 
S 


L 


2i6a  I  exp(-y2) 


a(y2  -  n20)  sin  gy  -  yVDt  cos(gy) 
a2(y2  ■  n20)2  +  y?Dt 


y 


where  0  °  A- 1  and  p  s  0a/  i/Dt.  This  now  permits  one  to  write 


f 

exp(v2  -  gv) 

dv 

J 

r(S,e) 

2  V  V2 

“  +  p  +  £lj 

V 

s 

r 

l 

2i0a  f  exp(-y2) 
^0 

|  a(y2  -  p,20)  sin( gy)  -  y-fDt  cos  gy 

2,  2 
a  (y 

-  p-20)2  +  y2Dt  J 

7 


Now  when  this  last  expression  is  put  into  the  expression  for  C  one 


gets  for  the  concentration  C 


00 


sin(p,Ti)  ~ 


OO 


Ii0a 


p.  sin(pr))  exp-p,  )dp. 


exp(-y2) 


a(y2  ~  M-2e)  sin(gy)  -  yV Dt  cos(gy) 
a2(y2  -  ii20)2  +  y2Dt 


3y 

y 


The  next  step  is  to  evaluate  the  first  integral  and  to  replace  p.  by  g, 
y  by  p.,  r)  by  y/  ^Dt  and  g  by  x/"V Dt  in  the  double  integral.  Then 


C  =  1  + 


liaB 


m 

C  00 

J  g  exp(-g2)  sin~==Jdg  J  exp(V 


a(n2  -  g20)  sin 

poc ' 

’ml 

-  p,VCt  cosj 

||XX  , 

mi 

l  H2 Dt  +  a2(jx2  - 

C2e)2 
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The  substitution  “  £/ ifi)t  and  jj.*  ■  p/VDi  gives  for  C,  after  removing  the 


primes j 


C  -l  +  ^x 


I  sin(£y)d£ 


«*  r  n  f  /  0 

f  ,  2  2s  a(u  " 

I  exp  -  U  +  H  )Dt  - ?p” 

n  L  ->  L  H  + 


a(u2  -  ?2e)  sin(px)  -  p,  cos(px) 


a  V  -  s2e)2 


.  (37) 


This  appears  to  be  Eq.  (1.9)  in  the  paper  by  Borisov  et  al.  ' 

The  following  formula  is  given  as  a  consequence  of  Eq.  (37)  with  some 
negligible  quantities  neglected. 


c  « i  -  ~  X 

Tt 

/r  2  sin(t)r)  +  cos(tf)  .  sl/2  1,2  ^  , 

I  - c~r - ^  sin  (zq)  '  u  exp(-t]  -  t)  z  e)  d  t)  a  z, 

-L  (1  +  z^)n  L  -I 


where 


VDt  ’  (Dt)1^2  Va6  * 


The  substitutions 


s20  .  ,n+  a/2 

z  =  a  -j2"  -  [i  and  q  =  p.(Dt) 


seer,  to  be  the  ones  used  to  get  Eq.  (38)* 

With  these  substitutions  the  integrand  in  Eq.  (37)  becomes 


Writ  .  hr  zp  2  rr  z-n  ^ 

,X  sin  VeDt  +  ea1^y  exp  _T1  "0  ■  0a 


2  Tj2  zrj  Vht  1  z  sin(rrj)  +  cos(rr)) 


(1  +  z  )ti 


31 


32 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  three,  major  published  works  on  grain-boundary  diffusion  have  been 
considered  and  the  following  conclusions  reached.  Whipple* s  solution  is  exact. 

The  solution  obtained  by  Borisov  et  al.  is  equivalent  to  that  obtained  by  Whipple. 
The  relationship  of  Fisher*  s  approximation  to  Whipple  *  s  solution  has  been  examined 
but  no  conclusion  about  its  validity  has  been  reached.  At  least  one  very  weak 
point  in  Fisher's  argument  has  been  found. 

Having  proceeded  this  far  one  may  recognize  rather  sharply  several 
specific  problems  that  could  very  profitably  be  considered.  Several  of  these 
are  enumerated  below. 

(1)  Numerical  Evaluation  of  Whipple* s  Solution.  With  the 
sure  knowledge  that  Whipple *s  solution  is  exact  it  seems 
appropriate  to  inquire  how  well  this  solution  describes 
experimental  results.  The  most  expedient  way  of 
accomplishing  this  would  seem  to  be  the  numerical 
evaluation  of  Whipple  *  s  solution.  Some  steps  in  this 
direction  were  taken  by  Whipple  himself,  however,  his 
work  does  not  cover  an  adequate  range  of  parameters. 

Furthermore,  his  calculations  are  based  on  a  formula 
which  results  from  taking  the  ratio  of  the  grain-boundary 
diffusion  coefficient  to  the  volume  diffusion  constant 

to  be  infinite  in  certain  of  the  places  where  it  occurs 
in  the  exact  solution.  This  approximation  is  probably 
quite  adequate,  however,  it  would  be  worthwhile  to 
evaluate  the  exact  solution  using  the  value  1CP  for  this 
ratio.  This  would  confirm  or  deny  our  present  speculation 
that  taking  &  *  a>  ,  as  Whipple  has  done,  does  not 
introduce  a  serious  error.  If  A  "»  in  a  satisfactory 
approximation  the  computational  effort  will  be  substantially 
reduced. 

(2)  Validity  of  the  Steepest  Decents  Approximation.  As  is 
often  the  case  with  steepest  decents  evaluations  of  integrals 
it  is  difficult  to  establish  the  range  of  validity  of  Whipple's 
approximate  formula.  One  possibility  is,  of  course,  the 
direct  comparison  of  the  approximation  with  the  exact  numerical 
results  discussed  above.  If  a  theorem  of  the  form  "If  the 
approximation  is  valid  for  a  concentration  c  then  it  is  valid 
for  any  concentration  less  than  c"  could  be  established 
analytically,  then  the  use  of  numerical  comparisons  to  determine 
the  range  of  validity  would  be  particularly  useful.  A  more 
ambitious  goal  and,  of  course,  one  which  is  less  likely  to  be 
realized  is  the  rigorous  analytical  determination  of  the  range 
of  validity  of  the  approximation. 

(3)  Analytical  Approximations.  A  rood  analytical  approximation 
to  Whipple's  solution  with  carefully  established  limits  of 
validity  would  be  extremely  valuable .  The  hope  of  finding 
such  an  approximation  is  rather  remote,  however,  the  search 
should  continue  as  a  part  of  any  sound  theoretical  program. 
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(1*)  Validity  of  Fisher 's  Approximation.  The  range  of 
parameters  for  which  Fisher’s  approximation  is  valid 
should  be  established.  This  investigation  should 
consider  the  following:  the  actual  deduction  of 
Fisher's  formula  using  Whipple’s  solution  to  study 
the  approximations  which  are  used;  the  shape  of  the 
isoconcentration  contours  compared  with  the  exact 
contours;  and  the  validity  of  the  activity  equation 
used  to  interpret  tracer  study  results. 

(5)  Borisov's  Deduction  of  Fisher’s  Formula.  Borisov^ 
claims  to  have  deduced  Fisher's  formula  from  the 
exact  solution.  His  presentation  is  sufficiently 
sketchy  that  it  requires  a  careful  reexamination. 

This  procedure  might  resolve  some  of  the  questions 
raised  above  and  at  any  rate  should  end  the  controversy 
which  has  been  going  on  for  several  years. 
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